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Abstract 

Since  departures  from  the  classical  assumptions  regarding  the  disturbances  in  a 
linear  regression  model  arise  frequently  in  empirical  applications,  several  computationally 
straightforward  procedures  are  presented  in  this  paper  for  testing  non— nested  models  when 
the  disturbances  of  these  models  follow  first—  or  higher— order  autoregressive  processes.  An 
empirical  example  is  used  to  illustrate  how  the  procedures  may  be  used  to  test  competing 
Keynesian  and  New  Classical  non— nested  models  of  unemployment  for  the  U.S.  using 
annual  time  series  data  for  1955-85. 


1.  Introduction 

Specification  tests  have  an  important  role  to  play  in  analyzing  the  adequacy  of 
econometric  models.  One  important  area  of  research  in  econometric  model  evaluation  is 
testing  non— nested  regression  models.  However,  research  on  this  topic  has  concentrated  on 
models  that  satisfy  the  classical  assumptions  of  serial  independence,  homoscedasticity  and 
normality  of  the  disturbances.  Since  departures  from  these  classical  assumptions  arise 
frequently  in  empirical  applications,  it  is  essential  to  develop  computationally 
straightforward  procedures  for  testing  general  forms  of  non— nested  econometric  models. 
By  a  general  form  of  a  model  is  meant  non-standard  cases  of  non— nested  linear  regression 
models  when  some  or  all  of  the  classical  assumptions  are  violated,  in  particular,  when  there 
is  serial  correlation  in  the  disturbances. 

The  paper  is  in  two  parts.  The  first  part  deals  with  problems  associated  with 
non— nested  models  when  the  disturbances  follow  first—  or  higher— order  autoregressive 
processes.  Since  residual  serial  correlation  may  arise  in  a  dynamic  specification,  lagged 
dependent  variables  are  permitted.  Several  asymptotically  valid  procedures  are  presented 
and  these  may  be  computed  from  standard  regression  packages.  An  empirical  example  is 
presented  in  the  second  part  of  the  paper  to  illustrate  the  use  of  the  non— nested  tests  when 
at  least  one  of  the  models  under  consideration  exhibits  autoregressive  disturbances. 

The  plan  of  the  paper  is  as  follows.  In  Section  2  we  discuss  alternative  asymptotic 
procedures  that  have  been  developed  in  the  literature  for  testing  non— nested  regression 
models,  and  examine  how  they  may  be  modified  to  take  account  of  alternative 
autoregressive  error  specifications.  The  relations  between  some  of  these  tests  are 
emphasized  to  simplify  the  steps  required  in  their  computation.  A  straightforward 
generalization  to  take  account  of  higher-order  autoregressive  error  specifications  is  then 
presented.  The  methods  developed  in  the  paper  are  then  used  in  Section  3  to  test 
Keynesian  and  New  Classical  models  of  unemployment  for  the  U.S.  using  annual  data  for 
the  period  1955—85.   Some  concluding  comments  are  given  in  Section  4. 
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2.   Test  Procedures  for  Models  with  Serial  Correlation 

In  the  case  of  non— nested  linear  regression  models  with  normal  and  spherical  errors, 
the  null  model  EL  is  tested  against  the  alternative  model  EL  and  the  two  models  are 
specified  as 

HQ:y  =   X/?+u0,       uQ   -    N(0,  a\  lQ) 

Hjiy  =  Z7+ux,      Uj   -  N(0,^In) 

in  which  y  is  the  n*l  vector  of  observations  on  the  dependent  variable,  X  and  Z  are 
n*k  and  n*g  matrices  of  observations  on  k  and  g  linearly  independent  regressors,  ft 
and  7  are  k*l  and  g*l  vectors  of  unknown  parameters,  and  u«  and  u,  are  n*l 
vectors  of  normally,  independently  and  identically  distributed  disturbances.  Strictly 
speaking,  normality  is  not  required  for  all  the  tests  to  be  discussed  below,  but  it  will  prove 

convenient  in  what  follows.    It  is  also  assumed  that  X  and  Z  are  not  orthogonal,  and  that 

— 1    '  — 1    '  — 1    ' 

the  limits  of   n     X  X,    n     Z  Z    and    n     X  Z    exist,  with  the  first  two  positive  definite 

and  the  third  non— zero.  If  X  and  Z  contain  stochastic  rather  than  fixed  elements,  the 
probability  limits  of  the  appropriate  matrices  must  exist,  and  X  and  Z  must  be  distributed 
independently  of  un   and  u, . 

In  order  to  take  account  of  departures  from  sphericality,  the  two  models  may  be 
rewritten  conveniently  as 


0:yt  =  xj/?+uQt  (1) 


H 

Hi  ;yt  =  z;-r+u1t 


i  •  't     "t ' r  ^t 

in  which  x'  and  z'  are  the  t'th  rows  of  X  and  Z,  respectively,  and  t  =  1,  2,  ...,  n.  In 
practice,  departures  from  the  classical  assumptions  can  arise  in  one  of  two  ways.  Economic 
theory  may  postulate  that  the  errors  in  a  model  follow  an  autoregressive  process  of 
specified  order.     Alternatively,  the  order  of  the  process  may  be  determined  empirically 


through  the  use  of  diagnostic  checks  and/or  information  criteria.  Regardless  of  the  ways  in 
which  the  presence  and  the  order  of  the  autoregressive  process  are  detected,  autoregressive 
errors  are  frequently  presented  as  an  integral  part  of  a  linear  regression  model.  Extending 
this  line  of  argument  to  the  non— nested  case,  such  models  may  themselves  display  different 
autoregressive  patterns. 

Pesaran  (1974)  derived  a  test  of  a  linear  regression  model  against  a  non— nested 
linear  alternative  when  the  regressors  of  both  models  are  non— stochastic  and  the  errors 
follow  stationary  first— order  autoregressive  (AR(1))  processes,  namely 

uit   =   Vit-l  +  £it  '     fit    -   NID(°.^)>     \Pi\<1  (3) 

for  i  =  0,  1  and  t  =  2,  3,  ...  n.  The  approach  was  based  on  a  direct  application  of  Cox's 
(1961,  1962)  procedure  for  testing  separate  families  of  hypotheses.  While  Pesaran's  test  is 
not  as  straightforward  computationally  as  some  of  the  procedures  to  be  discussed  below, 
and  modifications  are  required  to  account  for  lagged  dependent  variables,  it  does  have  an 
advantage  in  finite  samples  over  the  other  procedures  in  that  it  allows  for  stochastic  (rather 
than  fixed)  initial  values  of  the  AR(l)  processes  involved.  Of  course,  initial  conditions 
have  no  effect  on  the  properties  of  the  tests  asymptotically  and  would  be  of  concern 
primarily  in  finite  samples. 

Concentrating  solely  on  perceived  computational  advantages,  and  extending  the 
framework  of  straightforward  test  procedures  without  due  care  and  attention,  can  have 
some  undesirable  consequences.  A  simple  illustration  will  suffice.  There  have  been  several 
direct  applications  of  Davidson  and  MacKinnon's  (1981)  J  test  procedure  to  models 
exhibiting  serially  correlated  errors  (see  e.g.  Thornton  (1985)  and  Johannes  and  Nasseh 
(1985)),  although  the  original  test  was  derived  and  intended  for  models  with  serially 
uncorrected  errors.  A  J— type  test  of  H~  against  EL,  where  each  set  of  errors  obeys  (3), 
has  sometimes  been  interpreted  as  being  the  test  of  A  =  0  in  the  auxiliary  regression 
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yt   =  x;b  +  Ayu  +  vt  (4) 

in  which 

y«  =  z'o + ^iun-i  =  zv + 'i&t-i  -  zi-i?) 

is  a  consistent  estimate  of  the  predicted  value  of  y,  under  EL  and  b  =  (1  —  A)/?.  Since  v.  in 
(4)  corresponds  to  u^,  under  IL  :  A  =  0,  it  is  tempting  to  infer  that  a  direct  extension  of 
the  J  test  procedure  to  models  with  serial  correlation  is  valid,  namely,  presuming  that  v. 
follows  an  AR(1)  process  for  purposes  of  testing  EL  against  H,  via  a  test  of  A  =  0  in  (4). 

What  is  deceptive  about  such  an  approach?     Observe  that  EL  and  IL   may  be 
rewritten  as 

Ho :  h  =  Vt-i  +  (xt " Vt-i^  +  eot  =  *tCA  pq)  +  'ot  (5) 

Hi :  yt  =  ^yt-i  +  (zt  -^Vi)'7+  en  =  et(7,  px)  +  *lt  •  (6) 

A  linear  combination  of  (1)  and  (6),  with  weights  (1  —  A)  and  A,  respectively,  together  with 
the  addition  and  subtraction  of  Ay. . ,  leads  to  the  auxiliary  regression 

yt  =  (l-A)x;/?+Aylt  +  vt  (7) 

in  which 

vt  =  (l-A)uot  +  A[eu  +  (gt-ylt)] 

and  g.  =  gt(7,  Pi)-  It  can  easily  be  seen  that,  even  under  EL  :  A  =  0,  the  disturbance  term 
v.  in  (7)  will  be  asymptotically  correlated  with  y.,.  Observing  that  y,.  can  be  rewritten 
under  BL  as 

yit =  zi7+  Mxt-i^-zi-i^+ uot-i) 

and  denoting  the  probability  limit  of  p,  under  EL  by  plf)  (i.e.  plim^  p,   =  P-.0)>  it  is 

n-»tD 

straightforward  to  show  that 
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•    plimo  ln  \sn  yitvtl  =  Wio/(1_/,o)  • 

When  the  errors  under  EL  and  EL  are  serially  correlated,  y,.  and  v.  will  also  be 
correlated.  Thus,  if  serial  correlation  is  observed  in  competing  non— nested  regression 
models,  the  use  of  ordinary  least  squares  to  test  A  =  0  in  (7),  as  in  Backus  (1984)  and 
Milbourne  (1985),  would  lead  to  biased  test  statistics.  Moreover,  the  fact  that  y, ,  and  v. 
are  correlated  implies  that  an  efficient  method  of  estimating  p0  is  required.  As  a 
consequence,  the  Cochrane— Orcutt  procedure  will  not  be  valid  unless  it  uses  a  consistent 
estimator  of  p*  at  the  initial  stage.  An  appropriate  procedure  for  estimating  A  and  testing 
A  =  0  in  (4)  under  the  assumption  that  v.  follows  an  AR(1)  process  is  Hatanaka's  (1974) 
two— step  estimator,  which  is  valid  even  in  the  presence  of  lagged  dependent  variables  in 
both  models. 

An  alternative  procedure  for  testing  EL  against  EL  which  would  be  straightforward 
to  apply  in  the  present  context  is  to  test  EL  :  c  =  0  in  the  comprehensive  model  formed 
from  EL  and  EL  ,  namely 

yt  =x^b  + z^c  +  ut  (8) 

in  which  b  =  (1  —  A)/?,  c  =  A7,  and  u.  =  (1  —  A)ufi.  +  ^u,.,  which  follows  an  AR(1) 
process  under  IL  :  A  =  0.  It  is  now  well  known  for  the  case  where  both  ufi.  and  u,,  are 
serially  independent  that  the  standard  F  test  of  EL  :  c  =  0  in  (8)  yields  a  valid  test  of  EL 
against  H, .  Deaton  (1982),  Dastoor  (1983)  and  Gourieroux,  Monfort  and  Trognon  (1983) 
derived  a  non— nested  F  test  based  on  pseudo— true  values  of  selected  parameters  of  interest, 
McAleer  and  Pesaran  (1986)  showed  that  a  similar  analysis  could  be  conducted  using  Roy's 
union-intersection  principle,  and  Mizon  and  Richard  (1986)  based  their  F  test  on  the 
encompassing  principle.  However,  the  F  test  has  reduced  asymptotic  power  compared  with 
Cox— type  tests  (see  Pesaran  (1982a))  and  also  has  smaller  empirical  power  in  finite 
samples  (for  further  details,  see  Godfrey  and  Pesaran  (1983)).    The  generalization  of  the  F 
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test  of  c  =  0  in  the  case  of  serially  correlated  disturbances  is  a  much  more  complicated 
matter.  It  is  clear  that,  for  an  arbitrary  value  of  A,  the  composite  disturbance  u.  will  not 
follow  an  autoregressive  process,  and  maximum  likelihood  estimation  of  the  parameters  in 
(8)  will  not  be  a  straightforward  matter. 

On  the  basis  of  the  discussion  presented  above,  EL  may  be  tested  against  H,  after 
transforming  the  models  to  their  non— linear  counterparts  given  in  (5)  and  (6),  respectively. 
For  example,  it  is  valid  to  apply  the  comprehensive  model  approach,  or  the  procedures 
advanced  by  Davidson  and  MacKinnon  (1981)  and  Fisher  and  McAleer  (1981),  to  the 
non— linear  models  given  in  (5)  and  (6)  because  the  disturbances  are  serially  independent. 
Thus,  (7)  may  be  rewritten  as 

yt  =  (1  -  AJtyjyt,!  +  x{0-  pf^ft  +  Agt  +  e,  (9) 

where  e.  =  (1  —  A)c«t  +  Ae,..  The  unobserved  variable  g,  may  be  replaced  by  any 
consistent  estimator  of  its  components  under  EL.  An  obvious  candidate  is  g.  =  y, .  = 
z'  7  +  p-, u, ,_, ,  an  estimate  of  the  conditional  expectation  of  y,  under  EL ,  and  this  leads  to 
the  J  test  of  A  =  0.  If  y.^.  were  to  be  substituted  for  g,,  where  y,«.  denotes  the 
predictions  from  (6)  when  y,  is  replaced  by  yfi.  =  x'^+  PqU^.,  ,  the  predictions  from  BL, 
this  would  give  the  JA  test  of  Fisher  and  McAleer  (1981).  In  each  case  it  should  be  noted 
that,  unlike  (7),  the  disturbance  term  in  (9)  will  be  asymptotically  uncorrected  with  g. ,  or 
a  consistent  estimator  substituted  in  its  place.  However,  estimation  of  A  and  testing  A  =  0 
in  (9)  involves  non— linear  restrictions  between  the  parameters  (1  —  AW,  (1  —  A)/?, 
— pJl  —  \)j3  and  A.  This  means  that  a  test  of  A  =  0  in  (9)  generally  involves  non— linear 
estimation  and  the  testing  of  IL  by  this  method  on  standard  computer  packages  may  not 
always  be  possible. 

A  computationally  less  demanding  procedure  which  does  not  require  non— linear 
estimation  would  be  to  apply  the  P  test  of  Davidson  and  MacKinnon  (1981)  directly  to  (5) 
and  (6)  viewed  as  two  non— nested  non— tinear  regression  models.     The  appropriate  test 
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statistic  is  computed  as  the  usual  t-ratio  of  the  estimate  of  A  in  the  auxiliary  linear 
regression 

;ot = Fjd + A(;ot  -  ;lt) + <t  m 

in  which  £.,  =  y.  —  y..  denote  the  prediction  errors  under  H-  (i  =  0,  1),  F.  denotes  the 
maximum  likelihood  estimates,  under  HQ,  of  the  partial  derivatives  of  i^fi,  pQ)  in  (5)  with 
respect  to  0  and  pQi  namely 

Fj   =   (flft/»,  «t/»0)   =   [(xt-P0Xt-l)'^0t-l]  (11) 

U0t-1  =  yt-l~xt-l^ 
and 

d'  =  [(/3-^)',p0-p0]. 
Noting  that 

yt   =   y0t  +  e0t   =  ylt  +  6lt 

it  follows  that  efl.  —  c,.  =  y,f  —  y^,,  so  that  the  difference  in  the  prediction  errors  is  equal 
to  the  negative  of  the  difference  in  the  predictions  of  the  two  models.   Defining 

^10 1   E   6lt~  60t' 

substitution  of  (11)  into  (10)  yields,  after  some  algebraic  manipulation,  the  following 
convenient  auxiliary  linear  regression,  namely 

y*t  =  x^+<5uot-i-A7?iot  +  et  (12) 

in  which  y  =  y  -  PqYi  «,  x+.  =  x.  —  Pnxt_i  anc*  S  =  pQ  -  pQ.  The  expression  in  (12)  is 
didactically  more  appealing  than  that  in  (10)  because,  under  EL  :  A  =  0,  the  regression  of 
y^t    on    the    elements    of    x'      and    on    u«t _,     is    Hatanaka's    (1974)    two-step    (or 
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residual— adjusted  Aitken)  estimator  for  the  dynamic  adjustment  model  with  autoregressive 
errors.  In  the  context  of  this  problem,  then,  a  test  of  significance  on  77, «.  in  (12)  is 
equivalent  to  a  test  of  H^  given  in  (5)  against  EL  in  (6). 

A  similar  result  can  also  be  obtained  through  the  application  of  the  Cox  test  of 
Pesaran  and  Deaton  (1978)  directly  to  (5)  and  (6).   The  Cox  test  statistic  is  given  by 

I* 


No  -  VIW  (13) 


in  which 


TQ   =   (n/2)  log  (al/?10)  (14) 

V0(TQ)   =   (-o/-lo)^0-yio)'MF(yO-yio)]  <15) 

alo  =  ao  +  n-1(yo-yio)'(yo-yio)  (16) 

Mp  =   I-F(F'F)_1F' 

the  t'th  row  of  F  is  given  by  F'  in  (11),  and  y,  n  denotes  the  n*l  vector  of  estimates  of  the 
expectations  of  y,  under  H,  when  y,  is  replaced  by  yfi.  in  (6),  where  ynt  is  given  by 

yot  =  x^  +  ^u^  =  ^+Po(yt-rxt-ifl- 

Given  the  results  of  Lemma  3  of  Davidson  and  MacKinnon  (1981),  it  is  easily  established 
that  the  N0  test  in  (13)  and  the  tests  of  A  =  0  in  (10)  and  (12)  are  asymptoticall}' 
equivalent  under  BL. 

Since  it  is  asymptotically  valid  under  BL  to  replace  y,  n  in  (15)  and  (16)  with  y,,  an 
asymptotically  equivalent  expression  for  (14)  under  EL  is  given  by 


or 


TQ   =   (n/2)log[£jCl/(65e0+i7j:0»710)] 


TQ  =  -(n/2)  log  (1- 2c5J710/cie1)  (17) 
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in  which  17, «  =  e,  -  e«.   A  first-order  linearization  of  TL  in  (17)  is  given  by 

T0  =   efa/vl  .  (18) 

As  the  expressions  in  (14)  and  (18)  are  asymptotically  equivalent  under  HQ,  either  may  be 

substituted  into  (13)  to  yield  the  Cox  test  statistic.    However,  it  is  possible  to  derive  the 

1*  *  *  * 

asymptotic  distribution  of  n  * ^qViq  under  Hq,  and  to  base  a  Cox— type  test  on  cAi^n 

directly.   Suppose  (5)  and  (6)  are  rewritten  as 

HQ:yt   =  ft(AP0)  +  <0t   =  ft(")+£0t 
Hj  :  yt   =  gt(7,  Px)  +  «lt  =  gt(0  +  <lt 

with  maximum  likelihood  estimates  of  0'  =  (/?' ,  pfi)  and  0'  =  (7',  p,)  denoted  as  0  and 
(f>  under  EL  and  H,,  respectively.  If  (0' ,  <p')'  converges  in  probability  under  EL  to  (0x, 
0*)',  it  is  possible  to  rewrite  f^in  as 

^10  =  f6MF7?io +  V1' 

in  which  77,  ~  =  f(0«)  —  g(^+).  Under  suitable  regularity  conditions  (see  e.g.  White 
(1982)),  it  is  straightforward  to  show  that,  when  H~  holds, 

n    €6'io  ~  N(0,  aoplimn_l77ioMF77io)- 

n-^cD 

Therefore,  both  the  Cox  and  P  tests  are  asymptotically  equivalent  under  EL  to  the 
following  linearized  version  of  the  Cox  test,  namely 

LN0   =    ^o/I^ioM^lo)]*  (19) 

'2         —  1"   *  2 

in  which  ctq  =  n     eAeg  is  the  maximum  likelihood  estimate  of  cr~  under  EL  and  77,,,  = 

e,  —  6q  is  a  consistent  estimate  of  r/10  under  EL.    Notice  that  the  denominator  of  LNn  in 

(19)  can  be  computed  conveniently  as  the  residual  sum  of  squares  from  the  least  squares 

regression  of  77,  ^  on  the  columns  of  F. 
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The  two  tests  in  equations  (12)  and  (19)  would  be  expected  to  have  different 
properties  in  small  samples,  even  though  they  are  asymptotically  equivalent.  Recently, 
Bernanke  et  al.  (1988)  have  used  several  tests,  including  tests  which  were  developed  in  an 
earlier  version  of  the  present  paper,  to  test  non— nested  models  of  investment  subject  to 
serial  correlation.  On  the  basis  of  Monte  Carlo  experiments,  Bernanke  et  al.  (1988,  p.  320) 
suggest  that  the  P  test  given  in  equation  (12)  appears  to  be  the  best  of  those  presently 
available. 

The  results  obtained  above  can  be  generalized  straightforwardly  to  the  case  where 
the  disturbance  u.,  (i  =  0,  1)  follows  a  stationary  autoregressive  process  of  order  p.,  namely 
AR(Pi): 


Pi 
1=1  'iriH 


uit  "  . Viiuit-j  +  £ir     i  =  0-1  (2°) 


where  t  =  p+1,  p+2,  ...,  n  and  p  =  max(pn,  pA   In  this  general  setting,  the  P  test  of  EL 
is  simply  a  test  of  A  =  0  in  the  auxiliary  linear  regression 


y,t  =  x;^+jf/Juot-j-^iot  +  et 


in  which 


'0. 


*t  =  ^t-^/o/t-j 


Po 


*t 


"  xt-jVojVj 

'lOt   =   ylt  ~y0t 


-10- 


%   =   x^+^/0ju0t-j 
and 

.        Pi-     . 
Fit   =   ztTr+^/ljult-j' 

The  Cox  test  of  Pesaran  and  Deaton  (1978)  will  have  the  same  form  as  in  (13),  with  F^  in 
(11)  replaced  by 

Ft   =   Kt'  u0t— 1'  u0t-2'  -'  U0t-pQ)  • 

It  is  also  straightforward  to  obtain  a  linearized  Cox  test  as  in  (19)  corresponding  to  the 
case  where  u-,  follows  an  AR(p-)  process. 

So  far  in  this  section  we  have  considered  cases  in  which  the  disturbances  have  been 
serially  correlated.  The  situation  would  be  made  considerably  simpler  if  the  errors  were  to 
be  independent  but  not  identically  distributed.  Unless  the  actual  form  of  the 
heteroscedasticity  were  known,  transformation  of  (1)  and  (2)  along  the  lines  given  in  (5) 
and  (6)  (i.e.  in  the  case  of  known  autoregressive  structures)  would  not  be  possible. 
Nevertheless,  the  J,  JA  and  F  tests  of  EL  may  be  constructed  in  the  form  of  (4)  or  (8), 
where  the  form  of  heteroscedasticity  of  u.  is  unknown,  if  Wiute's  (1980) 
'heteroscedasticity— consistent  variance  estimator'  is  used.  Thus,  it  is  entirely 
straightforward  to  test  two  non— nested  linear  regression  models  against  each  other  in  the 
presence  of  general  forms  of  heteroscedasticity,  since  there  exist  computationally 
straightforward  non— nested  tests  which  require  estimation  only  of  an  auxiliary  linear 
regression. 

3.   Models  of  U.S.  Unemployment 

In  this  section  an  empirical  example  is  used  to  illustrate  the  application  of  the 
non-nested  P  and  LN  test  procedures,  as  given  in  equations  (12)  and  (19),  respectively. 
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The  example  centres  on  the  debate  between  the  New  Classical  model  of  unemployment  for 
the  U.S.  as  developed  in  Barro  (1977,  1979,  1981)  and  Rush  and  Waldo  (1988),  and  the 
Keynesian  (or  activist)  model  of  Pesaran  (1982b,  1988).  The  Keynesian  and  New  Classical 
models  are  given  as  follows: 

Keynesian  model:  Pesaran  (1988,  Appendix  Table  2) 

UNt  =  %  +  VxMILt  +  ^2UNt-l  +  ^3DMt  +  ^4DMt-l 

+  V>5DMt_2  +  V6t  +  V7WARt  +  etK  (21) 

New  Classical  model:  Barro  (1977),  Pesaran  (1982b,  1988),  Rush  and  Waldo  (1988) 
UNt  =  aQ  +  a1MILt  +  c^MINWj.  +  a3DMRHt  +  c^DMEH^ 

+  a5DMRHt_2  +  6tNC  (22) 

in  which  DMRH,  =  DM.  —  E._,(DM.)  is  the  error  term  in  the  money  supply  equation 
given  by 

DMt  =  PQ  +  ^1DMt_1  +  /52DMt_2  +  ^UN^ 

+  /?4Et_1(FEDVt)  +  DMRHt  (23) 

where  E^FEDV^  =  FEDVt  -  0.8DGRt  and  DGRt  =  DGt  -  E.^DG.)  is  the  error 
term  in  the  government  expenditure  equation  given  by1 

DGt  =  70  +  71DGt_1  +  72UNt_1  +  73WARt  +  DGRt.  (24) 

The  variables  in  equations  (21)  —  (24)  have  the  following  definitions: 
UNt  =  log[Ut/(l-Ut)] 
Ut  =  annual  average  unemployment  rate 
MIL,  =  measure  of  military  conscription 


1   Note  that,  in  defining  DGR,  by  (24),  it  is  implicity  assumed  that  the  value  of  WAR.  is 
known  to  economic  agents  at  time  t— 1,  that  is,  WAR.  is  perfectly  predictable  at  time  t— 1. 
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MINW.  =  minimum  wage  variable 

DM   =  rate  of  growth  of  money  supply  (Ml  definition) 

t 

DMRH    =  DM.  —  E,_,  (DM,)  =  unanticipated  rate  of  growth  of  money  supply 

I  LUX  I 

FEDV.  =  real  federal  government  expenditure  relative  to  its  normal  level 

E.    , (FEDV.)  =  anticipated  value  of  FEDVt  formed  at  time  t-1 

DG,  =  rate  of  growth  of  real  federal  government  expenditure 

DGR.   =  DG.  —  E.JDG.)  =  unanticipated  rate  of  growth  of  real  federal  government 

expenditure 
WAR,  =  a  step  dummy  variable  measuring  the  intensities  of  different  wars 
t  =  time  trend. 

The  primary  purpose  of  these  two  empirical  models  is  to  explain  the  rate  of 
unemployment  in  the  U.S..  Fiscal  and  monetary  variables  are  included  in  the  reduced  form 
Keynesian  model,  together  with  a  time  trend  to  explain  gradual  changes  in  the  natural  rate 
of  unemployment  over  time.  All  variables  in  the  Keynesian  model  are  observable  and  may 
be  efficiently  estimated  by  ordinary  least  squares  (OLS)  if  the  model  is  correctly  specified. 
However,  since  the  New  Classical  model  states  that  only  unanticipated  changes  in  the 
money  supply  affect  the  unemployment  rate,  a  sequential  estimation  and  testing  procedure 
may  be  used  to  generate  the  unobserved  money  supply  shocks  as  the  OLS  residuals  from 
the  money  supply  equation  in  (23).  Barro  (1977,  1979,  1981)  assumes  that  real  federal 
government  expenditure  relative  to  its  normal  level,  FEDV,,  can  be  anticipated  perfectly 
at  time  t-1  (that  is,  E,_,(FEDV  )  =  FEDV,).  When  this  unrealistic  assumption  is 
relaxed,  the  expected  value  of  FEDV  may  be  estimated  by  using  the  OLS  residuals  from 
the  government  expenditure  equation  in  (24).  The  variables  in  the  New  Classical  model 
include  monetary  shocks,  real  variables  to  explain  the  natural  rate  of  unemployment, 
countercyclical  monetary  and  fiscal  policy  variables,  and  variables  to  account  for 
government  financing  needs. 
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Maximum  likelihood  methods  may  be  used  to  estimate  the  New  Classical  model 
(that  is,  equations  (22)— (24)),  as  in  McAleer  and  McKenzie  (1989),  or  alternatively  a 
sequential  procedure  may  be  used,  as  in  virtually  all  of  the  papers  in  the  literature.  Using 
OLS  to  estimate  equations  (24),  (23)  and  (22)  in  a  sequential  manner  generally  yields 
inefficient  estimators  and  incorrect  standard  errors.  Pagan  (1984,  1986)  and  Pesaran 
(1987,  chapter  7)  provide  a  detailed  analysis  of  these  two  issues.  Murphy  and  Topel  (1985) 
also  examine  the  problem  of  incorrect  standard  errors,  and  McAleer  and  McKenzie  (1988) 
present  very  simple  proofs  of  several  of  the  existing  efficiency  results.  The  specific 
econometric  problems  associated  with  the  three— equation  New  Classical  model  are  analysed 
in  McAleer  and  McKenzie  (1989,  Appendix  A).  For  present  purposes,  it  is  sufficient  to 
state  that  sequential  OLS  estimation  yields  inefficient  estimators  and  the  standard  errors 
are  understated,  so  that  t— ratios  are  biased  towards  rejection  of  the  relevant  null 
hypotheses. 

Selection  of  the  sample  period  is  problematical.  Barro  (1977)  used  annual  data  for 
the  U.S.  for  1946-73,  Barro  (1979)  used  data  for  1946-77  and  1949-77  (the  latter  to  avoid 
problems  in  the  immediate  post— war  years  for  the  unemployment  equation),  and  Barro 
(1981)  chose  1946—78  (although  the  sample  period  became  1947—78  when  a  correction  was 
made  for  first— order  autoregressive  errors  in  the  unemployment  equation)  and  1949—78. 
Pesaran  (1982b)  used  data  for  1946-73;  Rush  and  Waldo  (1988)  chose  1946-73,  1949-73, 
1946—85  and  1949—85,  and  noted  some  evidence  of  serial  correlation  for  the  extended 
sample  period;  and  Pesaran  (1988)  used  1946—73  and  1946—85.  The  sample  period  used  in 
this  paper  is  1955—85.  This  covers  a  more  recent  period  than  that  considered  by  Barro  in 
his  three  papers,  and  also  avoids  the  major  dispute  between  Rush  and  Waldo  (1988)  and 
Pesaran  (1982b,  1988)  regarding  the  use  of  a  dummy  variable  for  war  by  the  former 
authors  to  accommodate  the  situation  whereby  the-  public  will  anticipate  and  has 
knowledge  of  the  quantitative  effect  of  an  abrupt  reduction  in  government  military 
spending  when  a  war  ends. 


-14- 


Both  estimation  and  testing  were  undertaken  using  the  computer  package  Microfit 
(see  Pesaran  and  Pesaran  (1989)).  The  OLS  estimates  of  the  Keynesian  model  are  given  in 
Table  1.  The  signs  and  magnitudes  of  the  estimated  coefficients  are  in  general  agreement 
with  those  of  Pesaran  (1988)  for  1946-85.  Since  the  sample  period  omits  the  effects  of 
World  War  II  and  the  Korean  War,  it  is  not  surprising  that  the  war  dummy  variable  is  not 
statistically  significant.  Moreover,  the  second  lag  of  the  money  supply  growth  rate  is  not 
significant.  Deleting  the  two  insignificant  variables  leaves  all  other  estimates  virtually 
unchanged  and  has  no  discernible  effect  on  the  outcome  of  the  non— nested  test  statistics 
reported  in  Table  3.  Finally,  there  appear  to  be  no  significant  departures  from  the  classical 
assumptions  of  correct  model  specification  or  from  normally,  independently  and  identically 
distributed  errors. 

Estimates  of  the  New  Classical  unemployment  equation  with  AR(1)  errors  are 
presented  in  Table  2.  The  government  expenditure  and  money  supply  growth  equations 
are  estimated  sequentially  by  OLS,  while  the  unemployment  equation  is  estimated  by  exact 
maximum  likelihood  subject  to  a  first— order  autoregressive  error  process.  Apart  from  the 
correction  for  serial  correlation,  the  estimated  magnitudes  and  signs  of  the  parameters  are 
very  similar  to  results  available  in  the  literature,  including  the  insignificance  of  the 
minimum  wage  variable.  The  diagnostic  tests  for  the  New  Classical  model  are  calculated 
using  the  adjusted  residuals  obtained  from  the  Cochrane— Orcutt  transformation  and  the 
first  derivatives  of  the  non— linear  function  resulting  from  the  AR(1)  errors.  The  classical 
assumptions  regarding  the  errors  and  correct  model  specification  appear  to  be  satisfied  for 
the  New  Classical  unemployment  equation. 

The  results  from  testing  the  non— nested  Keynesian  and  New  Classical  models 
against  each  other  using  the  P  and  LN  tests  are  reported  in  Table  3.  As  shown  in  Pesaran 
(1988)  and  McAleer  and  McKenzie  (1989),  the  non-rejection  of  the  Keynesian  model  is 
supported  strongly.  On  the  other  hand,  the  calculated  test  statistics  for  the  New  Classical 
model  are  beyond  conventional  critical  values,  leading  to  rejection,  although  it  should  be 
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noted  that  the  standard  errors  for  the  New  Classical  model  are  understated  so  that  the 
rejection  may  be  problematical. 

4.   Conclusion 

In  this  paper  we  have  presented  some  asymptotically  valid  and  computationally 
straightforward  procedures  for  testing  non— nested  regression  models  with  first-order 
autoregressive  disturbances.  The  procedures  were  also  generalized  to  take  account  of 
higher— order  autoregressive  processes  in  a  simple  manner.  An  empirical  example 
concerning  Keynesian  and  New  Classical  explanations  of  unemployment  for  the  U.S.  using 
annual  data  for  1955—85  was  presented  to  illustrate  the  use  of  the  testing  procedures. 
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TABLE  1 
OLS  Estimation  of  the  Keynesian  Model,  1955—85 


Regressor 

Coefficient 

Standard  error 

t— ratio 

INTERCEPT 

-2.0292 

.5018 

-4.0439 

MIL 

-4.0154 

1.1387 

-3.5264 

UN(-l) 

.3743 

.1255 

2.9815 

DM 

-5.4002 

1.7089 

-3.1600 

DM(-l) 

-3.5503 

1.8801 

-i.5477 

DM(-2) 

-1.7311 

2.2839 

-  .7580 

TREND 

.0365 

.0112 

3.2624 

WAR 

-.2123 

.2146 

-.9894 

R— squared 

.8945 

F— statistic  F(7,23) 

27.8537 

R— bar— squared 

.8624 

S.E.  of  regression 

.1134 

Residual  sum  of  squares 

.2958 

Mean  of  dependent  variable* 

-2.8340 

S.D.  of  dependent 

variable 

.3057 

Maximum  of  log— likelihood 

28.1187 

DW— statistic 

1.9735 

Durbins'  h— statistic' 

.1033 

Notes:        1   Thirty— one  observations  used  for  estimation  from  1955  to  1985. 
*   Dependent  variable  is  UN. 


Diagnostic  tests 

Test  statistics 

LM  version 

F  version 

(A)  Serial  correlation 

Chi-sq.(l)  =   .0043 

F(l,  22)  =  .0030 

(B)  Functional  form 

Chi-sq.(l)  =   .1063 

F(l,  22)  =  .0757  - 

(C)  Normality 

Chi-sq.(2)  =  2.2514 

Not  applicable 

(D)  Heteroscedasticity 

Chi-sq.(l)  =  1.9855 

F(l,  29)  =1.9845 

(A)  Lagrange  multiplier  test  of  residual  serial  correlation. 

(B)  Ramsey's  RESET  test  using  the  square  of  the  fitted  values. 

(C)  Based  on  a  test  of  skewness  and  kurtosis  of  residuals. 

(D)  Based  on  the  regression  of  squared  residuals  on  squared  fitted  values. 
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TABLE  2 
Estimation  of  the  New  Classical  Model  with  AR(1)  Errors,  1955-35 


Regressor 

Coefficient 

Standard  error 

t— ratio 

INTERCEPT 

-2.7629 

.2777 

-9.9501 

MIL 

-£.2247 

1.1653 

-5.3418 

MINW 

.3991 

.6849 

.5826 

DMRH 

-3.5209 

1.8827 

-1.8701 

DMRH(-l) 

-9.5880 

1.9896 

-4.8190 

DMRH(-2) 

-£.0384 

2.1697 

-2.7830 

R— squared 

.8509 

F— statistic  F(6,24) 

22.825: 

R— bar— squared 

.8136 

S.E.  of  regression 

.1320 

Residual  sum  of  squares 

.4180 

Mean  of  dependent 

variable* 

-2.8340 

S.D.  of  dependent 

variable 

.3057 

Maximum  of  log— likelihood 

22.6659 

DW— statistic 

1.9021 

Parameters  of  the  autoregressive  error  specification 


u  =      .4106   u(-l)  +  6 
(2.5073) 


t— ratio  based  on  asymptotic  standard  error  is  given  in  brackets 
Log-likelihood  ratio  test  of  AR(1)  relative  to  OLS:   Chi-sq.  (1)  =  5.1446 

Notes:       1   The  exact  inverse  interpolation  method  converged  after  5  iterations. 
2  Thirty— one  observations  used  for  estimation  from  1955  to  1985. 
*   Dependent  variable  is  UN. 
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TABLE  2  (continued) 
Diagnostic  tests 


Test  statistics 

LM  version 

F  version 

(A)  Serial  correlation 

Chi-sq.(l)  =   .0002 

F(l,  22)  =  .0002 

(B)  Functional  form 

Chi-sq.(l)  =   .5508 

F(l,  22)  =  .4115 

(C)  Normality 

Chi-sq.(2)  =  2.6699 

Not  applicable 

(D)  Heteroscedasticity 

Chi-sq.(l)  =  3.2371 

F(l,  28)  =3.3868 

(A)  Lagrange  multiplier  test  of  residual  serial  correlation. 

(B)  Ramsey's  RESET  test  using  the  square  of  the  fitted  values. 

(C)  Based  on  a  test  of  skewness  and  kurtosis  of  residuals. 

(D)  Based  on  the  regression  of  squared  residuals  on  squared  fitted  values. 


TABLE  3 

■ 

Non— nested  Test  Statistics 

- 

Null 
Model 

Alternative 
Model 

P  test 
Equation  (12) 

LN  test 
Equation  (19) 

Keynesian 

New  Classical 

-  .6948 

.7031 

New  Classical 

Keynesian 

3.5583 

-3.3197 

Note: 


The  P   and  LN   test   statistics  are  asymptotically  distributed  under  the  null 
hypothesis  as  N(0,  1). 
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